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J> ■ Abstract. We present a measurement of the mean intensity of the 

hydrogen-ionizing background radiation field at low redshift using 906 
Lya absorption lines in 151 quasar spectra from the archives of the 
£q ! Faint Object Spectrograph (FOS) on the Hubble Space Telescope (HST). 

■ Using a maximum likelihood technique and the best estimates possible 
H , for each QSO's Lyman limit flux and systemic redshift, we find J(vq)= 

7.6l|;o x 10" 23 ergs s" 1 cm" 2 Hz" 1 sr" 1 at 0.03 < z < 1.67. This is 

■ m good agreement with the mean intensity expected from models of the 
Q-f background which incorporate only the known quasar population. When 
q ■ the sample is divided into two subsamples, consisting of lines with z < 1 
H | and z > 1, the values of J(vq) found are 6.5lf 6 x 10" 23 ergs s 1 cm 2 

Hz" 1 sr" 1 , and l.Ojlo'l x 10~ 22 er S s s_1 cm" 2 Hz" 1 sr" 1 , respectively, 
indicating that the mean intensity of the background is evolving over the 
. ^ ' redshift range of this data set. Relaxing the assumption that the spec- 

^ . tral shapes of the sample spectra and the background are identical, the 

5-j | best fit HI photoionization rates are found to be 6.7 x 10" 13 s" 1 for all 

redshifts, and 1.9 x 10" 13 s" 1 and 1.3 x 10" 12 s" 1 for z < 1 and z > 1, 
respectively. 



1. Mean intensity of the UV Background 

The maximum likelihood method for measuring J{vq) as presented by Kulkarni 
& Fall (1993, KF93) consists of constructing a likelihood function of the form 




where 

f(N, z) = AN~P{1 + zy[l + cuiz)}-^. (2) 
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Figure 1. log[J(^o)] versus redshift: (filled triangle)- Shull et al. 
(1999); (filled squares & dotted line)- limits from Ha imaging; (crosses)- 
our results for z < 1 and z > 1; points with dotted error bars centered 
at z ~0.6,3.,4.5 are results from KF93, Paper II, and Williger et al. 
(1994), respectively; (solid curves)- HM96 models for two values of the 
global source spectral index, a 



Using the values of 7 and Aq from a separate maximum likelihood analysis, and 
a value of (3 from studies with high resolution data, eg. (3 = 1.46 from Hu et 
al. (1995), the search for the best-fit value of J(vq) consists of finding the value 
that maximizes this function, fixing the other parameters. 

The models of Haardt & Madau (1996, HM96) predict that the UV back- 
ground arising from QSOs drops by over an order of magnitude from z = 2.5 to 
z = 0. We therefore divide the sample into low and high redshift subsamples 
at z = 1 and solve for J{vq). These results, listed in Table 1. confirm some 
evolution in J(vq), though not at a high level of significance. The maximum 
likelihood analysis yields log[J(z/o)] = -22.18^q'q5 at z < 1 and log[J(i/o)] = 
-21.98^0 54 at z > 1. These results are also shown in Figure 1. 

Including associated absorbers, damped Lya absorbers, or blazars in the 
proximity effect analysis has little effect on the results. One might expect asso- 
ciated absorbers to reduce the magnitude of the observed proximity effect and 
hence cause J{vq) to be overestimated. The value found including the 45 asso- 
ciated absorbers in our sample is indeed larger, log[J(i'o)]=-21.74lQ '|g, versus 

log[J(^o)]=-22.1llo;40' but not significantly so. Likewise, if the intervening dust 
extinction in damped Lya absorbers is significant, including these objects in our 
analysis could cause us to overestimate the magnitude of the proximity effect 
and hence underestimate J{vq). However, the inclusion of these 7 objects only 
negligibly reduces the value of J{vq) derived. QSO variability on timescales less 
than ~10 5 years would be expected to smooth out the proximity effect distribu- 
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Table 1. Measurements of J{vq) 



Sample" Mines* lA 



P b log[J(i/ )] Qks 



1 259 0.8298, 6.73524 1 

la 162 1.5082, 4.92095 1 

lb 97 -0.8702,26.1886 1 

2 289 0.1502, 12.0134 1 




"Sample number- (1) All lines with W>0.32 A, (la) z<l, (lb) z>l; (2) All lines with W>0.24 

A 

^Number of Lya forest lines in sample 
C K-S probability 

tion (Bajtlik, Duncan, & Ostriker 1988). However, the inclusion of 6 blazars in 
the sample, all at z < 1, resulted in no discernible change in J(vq). The sample 
used in the analysis of HI ionization rates discussed below includes all of these 
objects. 

For each solution, we also execute a Kolmogorov-Smirnov (KS) test. The KS 
test provides a measure of how well the assumed parent distribution of lines with 
respect to redshift, given by Equation 2, reflects the true redshift distribution 
of lines. The KS probability, Qks> indicates the probability that a value of the 
KS statistic larger than the one calculated could have occurred by chance if the 
assumed parent is correct. The KS probability associated with each solution for 
J(vq) is listed in column 7 of Table 1. 

2. HI Ionization Rate 

Solving for the HI ionization rate, 



instead of for J(vq) avoids the assumption that the spectral indicies of the QSO 
and the background are identical. We modified our maximum likelihood code 
to conduct the search for this quantity and the results are listed in Table 2. 
Evolution in the UV background is more apparent in the HI ionization rate 
than in the solutions for J{vq). The result at z > 1 is 6.5 times larger than that 
at z < 1. 

We also parametrize the evolution of the HI ionization rate as a power law: 





T(z) = A(l + z) B 



(4) 



and solve for the parameters A and B in both the constant and variable threshold 
cases. The values we find are also listed in Table 2 and plotted in Figure 2. 
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Table 2. HI Ionization Rates 



Sample" 


j,A 


d 


b 


iogir m j 


1 


0.6925,7.64986 


1.46 


35 


i 9 -, 7 +U.5U 
-iz.i/_ 040 


la ... . 


0.8452,7.10998 


1.46 


35 


-12.70±°; 5 7 i 


lb ... . 


0.7209,7.29421 


1.46 


35 


-ii.88±g:« 


1 


0.6925,7.20637 


1.46 


35 


-12.67,1.7s 1 



a (l) All lines with W>0.32 A, (la) z<l, (lb) z>l 
1 Maximum Likelihood solution for A,B (see §2, Equ. 4) 
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Figure 2. HI ionization rate versus redshift: (points)- constant equiv- 
alent width threshold maximum likelihood solutions from this paper, 
z < 1 and z > 1 and from Paper II for 1.7 < z < 3.8; (dashed line)- So- 
lution to Equ. 4 for HST/FOS data alone; (solid line)- Solution to Equ. 
5 for HST/FOS data combined with high redshift data from Paper II, 
(dotted line)- HM96 solution to Equ. 5 
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HM96 parametrize their models of the HI ionization rate as a function of 
redshift: 

r(,)=^(l + ^exp^ " (z ^ c)2 ) (5) 

We combine our data set with that of Scott et al. 2000b (Paper II) to solve for the 
parameters A, B, zc, and S. We find (A, B, z c , 5)=(7.6 x 10~ 13 ,0.35, 2.07,1. 77), 
while the HM96 parameters for q = 0.5 are (6.7 x 10- 13 ,0.43,2.30,1.95). These 
results are shown in Figure 2, and the HM96 parametrization is also shown for 
comparison. 



3. Lya forest line density 

The number density evolution of Lya absorbers over the redshift range z = — 5 
cannot be approximated with a single power law. There is a significant break 
in the slope of the line number density with respect to redshift, near z ~ 1.7. 
Dave et al. (1999) show from hydro dynamical simulations of the low redshift 
Lyman a forest, that the evolution of the line density is sensitive mainly to the 
HI photoionization rate, but also to the evolution of structure. The flattening of 
dAf/dz observed by Weymann et al. (1998) is mostly attributed to a dramatic 
decline in T(z) with decreasing z. Dave et al. (1999) derive an expression for 
the density of Lyman a forest lines per unit redshift as a function of the HI 
photoionization rate: 

d -^ = C[{l + zfT^{z)]^H-\z\ (6) 

where C is the normalization at some fiducial redshift which we choose to be 
z = and T(z) can be expressed by Equ. 5. 

We fit the FOS and MMT absorption line data, binned in dAf/dz, presented 
in Dobrzycki et al. (2001, Paper IV) and in Scott et al. (2000a, Paper I), to 
this function in order to derive the parameters describing T(z) implied by the 
evolution in Lyman a forest line density. We observe some flattening of dAf/dz 
at z < 1.7, but not to the degree seen by Weymann et al. (1998) in the Key 
Project data. We find 7=0.50±0.21, for lines above a 0.24 A threshold (Paper 
IV), while Weymann et al. (1998) measure 7=0.15±0.23. We find (A, B, z c , S) 
= (3.0X10" 12 , 0.61, 5.5xl0" 7 , 7.07) and (1.9xl0 -11 , 0.38, 3.4xl0~ 7 , 6.21) for 
(Qm,^a)=(1.,0.) and lines with rest equivalent widths above 0.24 and 0.32 A, 
respectively. These fits to Equ. 6 are shown in Figure 3(a). In panel (b), we 
plot F(z), as expressed in Equ. 5, evaluated using the parameters found from 
the fit to Equ. 6 above. The HM96 solution and the solution derived from the 
full FOS and MMT data sets are represented by the thick and thin black lines 
respectively. The small values of z c derived from dAf/dz above translate into 
ionization rates that do not decrease dramatically with decreasing redshift and 
result from the less pronounced flattening of dAf/dz relative to the Key Project. 
The observed T(z) falls short of the ionization rate needed to fully account for the 
change in the Lyman a line density with redshift, indicating that the formation 
of structure in the low redshift universe plays a significant role in determining 
the character of the Lya forest line density. 
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Figure 3. (a) dM/dz versus redshift: H^ thr =0.24 A with fit to Equ. 
6 (solid points, dotted lines), IFthr =0.32 A with fit to Equ. 6 (open 
points, dashed lines), Equ. 6 evaluated with HM96 parameters for T(z) 
expressed by Equ. 5 (thick solid line), Equ. 6 evaluated with parameters 
for r(z) found in this paper (thin solid line); (b) T(z) versus redshift 
expressed by Equ. 5 using HM96 parameters (thick solid line), using 
parameters found in this paper (thin solid line), and using parameters 
found from fits to dM/dz for W thl =0.24 A and (O m ,^a)=(1.0,0.0) 
(thin dotted line), W thl =0.24 A and (fi M ,fi A )=(0.3,0.7) (thick dotted 
line), W thl =0.32 A and (n M ,n A )=(1.0,0.0) (thin dashed line), and 
W thT =0.32 A and (fi M ,n A )=(0.3,0.7) (thick dashed line) 



4. Summary of Results 

We have analyzed a set of 151 QSOs and 906 Lya absorption lines, the subset 
of the total data set presented in Paper III that is appropriate for the proximity 
effect. The primary results of this work are as follows: 

(1) The value of J(fo) is observed to increase with redshift over the redshift 
range of the sample data, 0.03 < z < 1.67. Dividing the sample at z = 1, we 
find J{uq)= 6.5lf g x 10 -23 ergs s" 1 cm~ 2 Hz -1 sr _1 , at low redshift and J(vq)= 
l.Oto'l x 10 -22 ergs s _1 cm~ 2 Hz -1 sr _1 at high redshift. 

(2) The inclusion of blazars at z < 1, damped Lya absorbers, or associated 
absorbers has no significant effect on the result. 

(3) Using information measured and gathered from the literature on each 
QSO's UV spectral index and solving for the HI ionization rate, yields 1.9 x 
10~ 13 s _1 for z < 1 and 1.3 x 10~ 12 s _1 for and z > 1. Solving directly for 
the parameters (A, B, z c , S) in the HM96 parametrization of T(z) using the 
HST/FOS data presented in Papers III and IV, combined with the high redshift, 
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ground-based data presented in Papers I and II, results in (A, B, z c , S)=(7.6 x 
1(T 13 ,0.35,2.07,1.77) for 1.7 < z < 3.8. 

(4) The z < 1 result is in agreement with the range of values of the mean 
intensity of the hydrogen-ionizing background allowed by a variety of local esti- 
mates, including Ha imaging and modeling of galaxy HI disk truncations (Mal- 
oney 1993, Corbelli & Salpeter 1993, Dove & Shull 1994, Kutyrev & Reynolds 
1989, Tumlinson et al. 1999). To within the uncertainty in the measurement, 
this result agrees with the one previous proximity effect measurement of the low 
redshift UV background (KF93). These results are consistent with calculated 
models based upon the integrated emission from QSOs alone (HM96) and with 
models which include both QSOs and starburst galaxies (Shull et al. 1999). The 
uncertainties do not make a distinction between these two models possible. 

(5) The results presented here tentatively confirm the IGM evolution sce- 
nario provided by large scale hydrodynamic simulations (Dave et al. 1999). This 
scenario, which is successful in describing many observed properties of the low 
redshift IGM, is dependent upon an evolving J{vq) which decreases from z = 2 
to z = 0. However, the low redshift UV background required to match the 
observations of the evolution of the Lyman a forest line density is larger than 
found from the data, indicating that structure formation is playing a role in this 
evolution as well. Our results and the work of others are summarized in Figure 
1. We find some evidence of evolution in J(fo), though it appears that even 
larger data sets, especially at z < 1 and/or improved proximity effect ionization 
models will be required to improve the significance. 
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